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The durations and average speeds of ultrashort optical pulses transmitted through chiral sculptured thin films (STFs) were calculated 
using a finite— difference time-domain algorithm. Chiral STFs are a class of nanoengineered materials whose microstructure comprises 
parallel helicoidal nanowires grown normal to a substrate. The nanowires are ~ 10-300 nm in diameter and ~1-10 /xm in length. Durations 
of transmitted pulses tend to increase with decreasing (free-space) wavelength of the carrier plane wave, while average speeds tend to 
i: increase with increasing wavelength. An increase in nonlinearity, as manifested by an intensity-dependent refractive index in the frequency 
■ domain, tends to increase durations of transmitted pulses and decrease average speeds. The circular Bragg phenomenon exhibited by 
' a chiral STFs manifests itself in the frequency domain as high reflectivity for normally incident carrier plane waves whose circular 
polarization state is matched to the structural handedness of the film and whose wavelength falls in a range known as the Bragg regime; 
films of the opposite structural handedness reflect such plane waves little. This effect tends to distort the shapes of transmitted pulses 
with respect to the incident pulses, and such shaping can cause sharp changes in some measures of average speed with respect to carrier 
, wavelength. A local maximum in the variation of one measure of the pulse duration with respect to wavelength is noted and attributed 
. to the circular Bragg phenomenon. Several of these effects are explained via frequency-domain arguments. The presented results serve 
i: as a foundation for future theoretical and experimental studies of optical pulse propagation through causal, nonlinear, nonhomogeneous, 
and anisotropic materials. 



1 Introduction 

An ultrashort optical pulse is a pulse of light of duration less than ~10 optical periods, where the optical 
period is calculated from the dominant wavelength in the pulse spectrum. Such ultrashort pulses are 
generated with commonly implemented mode-locked lasers, and several schemes for compressing and 
amplifying them have been formulated [1]. Optical pulses are used for micromachining [2], lidar [3,4], and 
quantum control and probing of wavepackets and chemical reactions [5,6]; and their use has been suggested 
for telecommunications [7]. 

It is necessary to understand how ultrashort optical pulses propagate in materials in order to design 
new materials and devices to manipulate them. Earlier, we reported on the spatiotemporal evolution of 
ultrashort optical pulses as they propagate through complex nano-engineered materials called sculptured 
thin films (STFs). A STF is a unidirectionally nonhomogeneous, anisotropic medium whose morphology, 
which consists of an array of parallel shaped nanowires, is engineered by design [8]. The morphology 
of STFs may break certain translational and rotational symmetries, and these films possess dispersive 
dielectric and possibly dispersive magnetic properties in the optical regime. 

The circular Bragg phenomenon is exhibited by chiral STFs due to their periodic nonhomogeneity and 
structural chirality arising from the constituent nanowires being nanohelixes as shown in the electron 
micrograph of Figure ^ As circularly polarized monochromatic light propagates through vacuum, the 
electric field traces out a helix. If the handedness of this helix matches that of the nanohelixes in a chiral 
STF, and if the wavelength of the light lies within a range of wavelengths called the Bragg regime, the 
light normally incident on the chiral STF is mostly reflected; if the handedness of the helix is opposite, the 
reflection is little. In the time domain, the circular Bragg phenomenon is manifested as a pulse-bleeding 
phenomenon [9,10]. 

Previous research has shown that chiral STFs can affect the shapes, durations, and energy content 
of ultrashort optical pulses depending on the circular polarization state of the incident carrier plane 
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Figure 1. Scanning electron micrograph of a chiral sculptured thin film; note the periodic helical morphology (Courtesy: M. W. Horn, 

The Pennsylvania State University). 

wave [10]. As the thickness of a chiral STF increases, the circular Bragg phenomenon ripens: the reflection 
first increases and then saturates, and differences in shape and duration between incident and transmitted 
pulses arise [11]. Ultrawideband pulsed plane waves, known as videopulses, ultranarrowband rectangular 
pulsed plane waves, and two-dimensional pulsed beams can therefore also be shaped by chiral STFs [12-14]. 
Pre- and post-resonance Bragg regimes have also been elucidated [15]. 

So far, all research into the shaping of pulsed plane waves and pulsed beams by chiral STFs has been of 
a qualitative nature. The shapes of pulses reflected, refracted, transmitted, and diffracted by chiral STFs 
were predicted numerically, but measures of the extent of pulse shaping were not introduced. Our purpose 
here is to apply three measures of the pulse duration — the equivalent, root mean square, and correlation 
durations — and three measures of the average speed of the pulse — the peak speed, the center-of-energy 
speed, and the correlation speed — to the shaping of pulsed plane waves of different total energies, carrier 
wavelengths, and carrier polarization states after transmission through linear and nonlinear chiral STFs. 

This work is important for four reasons. First, the measures of pulse shaping will allow direct comparisons 
between our theoretical conclusions and future experiments. Second, quantification of pulse shaping by 
complex materials like STFs will facilitate the design of STFs to shape pulsed plane waves and pulsed 
beams in controlled ways. Third, the algorithms we develop can be used to evaluate designs of STF-based 
devices; and fourth, they advance the state of the art of the time-domain analysis of electromagnetic pulse 
propagation to include propagation in causal, inhomogeneous, anisotropic, and nonlinear materials. 

Parenthetically, by virtue of mathematical isomorphism we expect our results to hold, in general, for 
cholesteric liquid crystals (CLCs) too. We note that others have modeled optical propagation in noncausal 
nonlinear cholesteric liquid crystals (CLCs) in the frequency domain [16] and with the slowly-varying 
envelope approximation [17]; propagation in linear CLCs has also been undertaken with the finite element 
method [18]. Second harmonic generation in nonlinear chiral STFs was also studied previously in the 
frequency domain [19]. Our work distinguishes itself from these earlier studies on account of its causal 
time-domain treatment of the propagation of pulsed plane waves through nonlinear chiral STFs without 
the invocation of the frequency domain or slowly varying envelope approximation. 

Our report on the quantification of the shaping of pulsed plane waves by chiral STFs proceeds as 
follows. The constitutive relations of a chiral STF, along with a finite-difference scheme for calculating the 
propagation of pulsed plane waves through vacuum and chiral STFs, are presented in Section [21 Measures 
of the durations and average speeds of pulsed plane waves are also introduced in Section EJ We discuss 
the main results of our work in Sectional and present the main scientific and technological conclusions in 
Section ^ 



2 Theory 

This section describes the constitutive relations of spatially local, dielectric, structurally chiral STFs; our 
choice of ultrashort optical pulses; the numerical method we use to predict their propagation; and three 
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measures each of pulse duration and speed. 

2.1 Maxwell curl postulates and constitutive relations 

We begin with the source-free time-domain Maxwell curl postulates 

V X E{r,t) = -dtB{r,t) , (1) 

V X H{r,t) = dtD{r,t) , (2) 

where and B represent the primitive macroscopic electric and magnetic fields, respectively, while D and 
H are the corresponding induction fields. The electromagnetic field is a function of r, the position vector, 
and t, the time. We employ a cartesian coordinate system defined by the set of unit vectors {u^,Uy,Uz}- 
The time derivative is denoted by dt- 
A spatially local, dielectric, structurally chiral STF possesses the constitutive relations 

mr,t) = eoE{r,t) + P{r,t) , (3) 
B{r,t) = ,ioH{r,t), (4) 

where eo is the permittivity and //q is the permeability of vacuum. The polarization P is related to the 
electric field by the following equation: 



P {r, t) = eo r (r, t') ■ (l + p^i \E (r, t-t')f)E (r, t - t') dt' . 

J 



(5) 



When pni = 0, the film exhibits a linear response to an applied field. When pni ^ 0, the response is 
third-order nonlinear; specifically, the material exhibits an intensity-dependent refractive index [19,20]. 
The nonlinearity is thus isotropic. Whether linear or nonlinear, the chosen material does not respond 
instantaneously to the applied field. 

The dielectric susceptibility dyadic x (zii t) contains the information about the nonhomogeneous, 

anisotropic, and temporally nonlocal properties of a chiral STF. Suppose the chiral STF occupies the 
region [zl,zr]. Then the dielectric susceptibility dyadic, which is null-valued in the region z ^ [zL,zji], 
can be factored as: 

^ (l, t) =g{z- Zl) ■ ^^^^ {t) ■ {z-Zl) , ze [ZL, Zr] , (6) 

where 5 is a rotation dyadic that describes the morphology of the chiral STF along the z axis (the axis of 
nonhomogeneity) , and x (^) is a reference susceptibility dyadic that describes, in accord with causality, 

=ref 

the time-dependent and anisotropic response of the material. In this paper we treat unidirectionally nonho- 
mogeneous STFs, and therefore exclude the recently fabricated bi- and tri- directionally nonhomogeneous 
films from consideration [21,22]. 
The rotation dyadic takes the form 

g{z)=g^{z)-gy{a) , (7) 

where the dyadic 

(^) = cos (J^^ [u^u^ + UyUy) + /isin [uyU^ - u^Uy) + u^u^ (8) 
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depends on the structural half period and the structural handedness parameter /i, and the dyadic 

(a) = cos a {u^Ux + ILzl^z) + sin o {Uzlhc ~ HxlLz) + iLylLy (9) 

depends on the angle of rise a. The chiral STF is either structually right-handed when /i = 1 or structurally 
left-handed when h = —1. 

The reference susceptibility dyadic, which takes on an orthorhombic form due to the physical vapor 
deposition process by which chiral STFs are made [8], is defined as 

X Ai) = Xi (t) u^Ux + X2 {t) u u + X3 {i) UzUz , (10) 

=rej " " 

where 

Xb (*) =Pfi'^Bexp sin(LJ6t)Z^(t) , i3 = l,2,3, (11) 

is determined both by the oscillator strengths and by the parameters Ajj and A'jj, which quantify the 
resonance wavelengths and absorption characteristics of the film; here, cjb = 27rco/AB are the resonance 
angular frequencies and U (t) is the unit step function. The Lorentzian dependences Hll|) are based on 
a classical material model that treats the charges in a material as undergoing damped oscillations when 
disturbed from equilibrium by an electric field [23]. 



2.2 Matrix partial differential equation 

We consider pulsed plane waves incident normally on chiral STFs, and hence there is no variation in the 
electromagnetic field along the x and y directions — i.e., dx = dy = 0. Upon writing the components of 
the electromagnetic field in a 6-vector as [F_{z,t)] = [E^, Ey, Ez, H^, Hy, Hz]^ , where the superscript 
indicates the transpose, and substituting the constitutive relations into the Maxwell curl postulates © 
and 0, we obtain the matrix partial differential equation 



dt [F {z, t)] =co\Y^dz [F {z, t)] - (1/eo) dt \Q (z, t) 



(12) 



In this equation, 



m 



-r?o 

% 



1/r/o 

-1/r/o 





(13) 



is the vacuum propagation matrix, r/o = y^/xo/eo being the impedance of vacuum and cq = l/^eo/^o the 
speed of light in vacuum, and the column vector 



[Q {z, t)] = [P, {z, t) , Py {z, t) , Pz {z, t) , 0, 0, 0]'^ 



eo 



1^ (l+Pni\E{z,t-t')f) [m{z,t')] [F{z,t-t')]dt' . 



(14) 
(15) 
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In the previous equation, 



[m(z,t)] 



=e J — 



(16) 



contains the linear material properties, with [£] as the 3x3 null matrix. We denote the matrix equivalents 
of vectors and dyadics by enclosing the corresponding symbols in square brackets. 

Equation (|12|) can only be solved approximately, analytical solutions being virtually impossible for the 
materials being considered here. It was discretized with respect to both z and t, using the leapfrog method. 
The details are available elsewhere [24]. 



2.3 Incident pulsed plane waves 

We chose the incident signal as a carrier plane wave that is amplitude-modulated by a pulse envelope 
ip (t), i.e., at z = we prescribe the boundary condition 



[F(0,t)]=V'(t) 



cos (uJcart) 
COS (uJcart + 4>) 



(1/r/o) cos {Ucart + 
(1/%) cos {iVcart) 





U{t) 



(17) 



where ujcar = '^T^c^/Xcar IS the carrier angular frequency and Xcar is the carrier wavelength, while (p 
is a carrier phase. When (p = 7r/2, the carrier plane wave is right circularly polarized (RCP); when 
(j) = — 7r/2, the carrier plane wave is left circularly polarized (LCP). The initial condition was chosen to be 

[Z(^,o)] = [o]. 

Even though we treat pulsed plane waves in this paper, it should be noted that no electromagnetic 
signal — or STF — can have an infinite transverse extent; indeed, all physical electromagnetic signals are 
pulsed beams. However, the pulsed-plane-wave approximation to a pulsed beam is appropriate when the 
beamwidth is much larger than the carrier wavelength (but still smaller than the transverse extent of the 
STF). The characteristics of two-dimensional pulsed beam propagation in chiral STFs have been reported 
elsewhere [14]. 

We chose the pulse envelope ip (t) to be gaussian, so that 



iP{t) 



VoUi 



exp 




t>0, 



(18) 



but tp (t) = for t < 0. The time constant tq governs the duration of the pulse, while is the time delay 
which permits the whole pulsed plane wave to enter the computation domain. The prefactor on the right 
side of H18() is chosen so that 



S,{z,t) dt = Ut 



(19) 



Sz = Ikz' {E.^ H_) being the axial component of the instantaneous Poynting vector. Hence, Ut is the total 
energy per unit area of the incident pulsed plane wave. 
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2.4 Quantification of pulse durations and average speeds 

We now turn to quantitative measures of pulse duration and average speed. Our goal is to define measures 
that are useful for experimental research or theoretical insight or both. As the chosen pulses are highly 
broadband and suffer appreciable distortion, the usual concept of group velocity [25] is not meaningful for 
analysis except in frequency-domain arguments. More sophisticated measures were therefore devised. 

The pulse duration is an important parameter because it determines the time scale for measurements 
or manipulations made with the pulse; pulses of shorter duration permit the measurement of events of 
shorter duration [26]. In order to define measures of pulse duration, we considered the m}^ moment -M^f^ 
of a scalar function / (^) with respect to variable ^ on the interval [^a, 



A^^^^ (/ (0 , [u ^b]) = J^^ r / io ■ (20) 

The root mean square (RMS) deviation af^^ from the centroid 
of / (0 "^ith respect to ( on the interval [^ajCfe] is defined as follows [27]: 



c/,€ (/ (e) , [Ca, ^b]) = (21) 



f\^-Cf,d'fiOd^, (22) 



A^g(/(6,[e.,6]) 

Both and ct/^^ are undefined when if (0 ^ i^a, Cb]) = 0- 

One measure of the pulse duration, as evaluated over a time interval [ta,tb\ at a point Zr, is 

Tu = 2au,t{U{Zr,t),[ta,tb\) , (24) 

where the electromagnetic energy density 

U iz,t) = eomz,t)\'' + no\H{z,t)f . (25) 

Alternatively, we could compute the RMS deviations from the centroid of 5^ {z, t) or the intensity \E_ (z, t)f 
with respect to time. Another measure of the pulse duration, called the equivalent duration [27], is ■M^\ 
computed at a particular Zr over an interval [taj^b]) and divided by the peak value of U{zr,t) over that 
interval, i.e., 

^ M^^,{U{Zr,t),[taM) 

Msx{U{Zr,t),[ta,h]) ' ^ ^ 

where Max indicates the maximum value of a function over the chosen time interval. The correlation 
duration Tc is defined by [27] 

IZoIZoUi^r,t)UiZr,t + t')dt'dt 

"^c = M ■ (27) 

MP.,{UHZr,t),i-^,^)) 
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We decided to use Tp, Tu, and Tc as measures of pulse duration in this study. 

There are many measures of the pulse speed, including the group speed and the centrovelocity [28]. We 
calculated three measures of the speed of the pulse through the film: the average speed of the pulse's peak 
value of U {zr,t), the average speed of Cu,t, and the average correlation velocity [29]. The energy density 
U {zr, t) of the transmitted pulse was recorded over the interval [ta, tb] at some location Zr > zr. Then the 
time 

tp = tMdJKiUiZr,t),[ta,tb]) (28) 

where the function ty^g^ indicates the time at which the maximum value of its argument occurs, the time 

tu = Cu,t{U{Zr,t),[ta,tb]) , (29) 

and the time 



tc = td + U (0, t') U{zr,t + 1') dt'^ 



(30) 



were found. From these data, the average speeds Cp, Cu, and Cc were calculated from the formula 

CB Zr - Zl 



CO Co (tfi - td) - {Zr - (zR - zl)) ' 



B = p,u,c. (31) 



Thus the speeds Cp, c„, and Cc are the thickness of the chiral STF divided by the computed times tp, tu, 
and tc, less the time over which the pulse traverses the vacuous regions < z < z^ and zr < z < Zr- We 
note that these measures are average speeds because they are in the form of a finite distance over a finite 
time interval, and hence do not reveal variations in the speed of the pulse as it enters, traverses, and exits 
the chiral STF. For the foregoing calculations, wc approximated integrals with the rectangular rule [30]. 

The measures of pulse duration and average speed have different physical meanings and concomitant 
advantages and disadvantages. The equivalent duration Tp is the duration a rectangular pulse would have, 
if it possessed the same peak value of U {zr,t) and the same total energy per unit area as the pulse 
under consideration. Its advantage lies with its simplicity; however, it depends crucially on the correct 
determination of the peak value of U {zr,t). This same problem is a disadvantage of Cp, which is the 
average speed of the peak of the pulse through the chiral STF. This measure of the speed can even become 
multi-valued if the peak value of U {zr,t) occurs at more than one local maximum (like the twin humps 
of a Bactrian camel). 

The duration is the RMS deviation of U {zr, t) from its first moment, normalized to the total energy 
per unit area of the pulse. It is analogous to the concept of radius of gyration in mechanics. We call c„, the 
average speed of C^.i, the center-of-encrgy speed. It is analogous to the concept of the velocity of a body's 
center of mass in mechanics, and is closely related to the centrovelocity and a redefinition of the group 
velocity^, as discussed elsewhere [28]. The advantage of these measures is that they depend on U{zr,t) 
over the whole pulse, and not just on the peak value of U {zr,t). However, it is possible for the center of 
energy to be located at a value of z at which U {z, t) = [29] . 

The correlation duration Tc is a measure of how tightly U {zr,t) is concentrated in time. The correlation 
velocity Cc is the average speed of the peak of the correlation of the incident pulse with the transmitted 
pulse. Like Cp, it depends heavily on the shape of the transmitted pulse. In some cases, it could become a 
multi-valued function of Zr and the chosen time interval [ta,tb]- However, it does provide a way to track 
the propagation of energy in the pulse in a different way than c„, and it is fairly robust [29]. 

All six of the quantities we have defined in this subsection are amenable to experimental determination, 
for example by frequency-resolved optical gating (FROG) [26]. In our definitions of the measures, we 

^The centrovelocity is the velocity of C|b|2 t (jKl'^ (zr,t) , [tajtt]^ [28]. 
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have chosen to use U {zr,t) instead of \E{zr,t)\'^, for example, because the former is a Lorentz-invariant 
quantity. 

3 Results and discussion 

We now turn to a discussion of the results of our calculations. 



3.1 Constitutive and excitation parameters 

The chiral STF we studied possesses the following constitutive parameters: pi = 0.52, p2 = 0.42, = 0.40, 
AT-^ 2,3 = 100, Ai = 290 nm, A2,3 = 280 nm, O = 200 nm, and a = 20°. The film is structurally right- 
handed, i.e., h = +1, and we chose it to be either linear (i.e., Pni = 0) or nonlinear (pni = 3 x 10~^^ m^/V^), 
in accord with simple order ~of-magnitude estimates [31]. We fixed zl = 30 /xm and Zfi = 34 /im; thus, the 
chiral STF is ten pitches thick. The variable z was discretized into intervals of Az = 2 nm, and the variable 
t into intervals of At = 6.34 x 10~^ fs. The center wavelength of the Bragg regime of the linear chiral STF 
is approximately 516 nm, and its full-width half-maximum bandwidth is approximately 27 nm [10,32], 
though we note that both these quantities may change slightly when the chiral STF is nonlinear [33]. For 
all the results presented here, the incident pulse has the time constant tq = 2 fs and time delay trf = 8 fs. 
We calculated for three different combinations of material nonlinearity and incident pulse energy per unit 
area: 

Case L: a linear film, with Ut = 1 x 10^ J/m^, 
Case Nl: a nonlinear film, with Ut = 1 x 10^ J/m^, or 
Case N2: a nonlinear film, with Ut = 2 x 10^ J/m^. 

Two indexes of refraction that govern the speed of monochromatic plane waves through a linear chiral 
STF are [8] n2 (Aq) = y^i2 (Ao) and (Aq) = y^e^ (Aq) where 



€i (Ao) = l + xe (Ao) = 1 + 







1 + (27riV,)-2 


( 2-KCa ' 








1 + (27riV,)-2 


( 27rCo \ ^ 


2 2 

1 ( 2cQUJt\ 



1,2,3, (32) 



and 



erf (Ao) 



h (Ao)e3(Ao) 



ei (Ao) sin^ a + es (Ao) cos^ a ' 



(33) 



where " indicates the temporal Fourier transform and Ao is the free-space wavelength. Plots of the real and 
imaginary parts of n2 (Ao) and (Ao) are shown in Figure |21 The corresponding phase velocities 

Vp2 = Co/Re (722) , (34) 
Vpd = Co/Re [rid) , (35) 



and group velocities 



Vg2 = Co/Re (n2 + ujd^n2) , (36) 
Vgd = Co/Re {nd + ujOu^rid) , (37) 

are shown in Figure |31 here, a; = 27rco/Ao and d^^ = d/dto. These frequency-domain plots help to explain 
our time-domain results. 
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3.2 Reflected and transmitted pulses 

The spatiotemporal evolution of the incoming pulse was computed on the spatial domain z S [0, 60] fim, 
of which the chiral STF occupies z £ [30, 34] fim. Measures of pulse duration and speed were computed 
at Zr = 36 /im over the time interval [0, 190] fs. The wavelength Xcar of the carrier plane wave was set to 
either 415, 515, or 615 nm; and the polarization state of the carrier plane wave was chosen as either LCP 
or RCP. At visible wavelengths, the selected linear constitutive properties yield refractive indexes that are 
commonplace in optics [34]. 

Let us now examine the effects of the chiral STFs on pulses reflected from and transmitted through 
them. We are interested mainly in trends in the pulse shape, duration, and speed as the pulse energy per 
unit area, carrier wavelength and polarization state, and the nonlinearity of the chiral STF are varied, 
as opposed to particular values that the duration, etc., may take on as those parameters are varied. 
The shapes of the incident pulses, as exemplified by snapshots of their energy density in Figure [IJ are 
independent both of the carrier wavelength and carrier polarization state. The incident pulses traversed 
vacuum (z G [0,30] /um), were scattered by the STF, and eventually became transmitted signals in the 
vacuous region z G [34, 60] /im and reflected signals in the region z G [0, 30] fim. 

Let us begin with Case L. Both transmitted and reflected pulses are captured in the snapshots of U {z, t) 
at time t = 190 fs in Figure |S1 where an important trend can be perceived. Irrespective of the polarization 
state of the carrier plane wave, the peak energy density of the transmitted pulse increases with increasing 
carrier wavelength. That is, the peak value of U {z,t), z > zr, is greater when Xcar = 615 nm than when 
Acar = 515 or 415 nm for a given carrier polarization state. This trend applies in reverse to the lengths (i.e., 
z-extents) of the transmitted pulses: the transmitted pulses are longer when Xcar takes on lower values. 

The effects of structural handedness are most evident in the plots of the reflected pulses shown in 
Figure ini which contains plots of U {z, t) as a function of time at z = 27.6 /im. The incident pulse is either 
LCP or RCP, and Xcar = 415, 515, or 615 nm. The energy densities of the reflected pulses (which begin at 
t ~ 110 fs) are much smaller than those of the incident pulses. 

The effects of the circular Bragg phenomenon are evident in Figure El In the bottom row of plots, for 
which the conditions for the circular Bragg phenomenon are satisfied [10], the primary reflected pulse 
possesses both a larger head and longer tail than when the incident carrier plane wave is LCP (top row 
of plots). The longer tail is due to the bleeding of energy from the refracted pulse inside the film to the 
refiected pulse [10], and it is most pronounced when Xcar = 515 nm, i.e. when the carrier wavelength lies 
in the Bragg regime. 

A secondary pulse appears in all six cases beginning at t ~ 155 fs. This smaller pulse is the result of the 
refiection from the chiral STF / vacuum interface at z = zr. When the incident light is RCP, U {z,t) at 
z < zl oscillates rapidly in both the head of the primary pulse and the secondary pulse, a phenomenon 
that is due to interference between LCP and RCP components of the reflected light [33]. The tail of the 
primary reflected pulse does not exhibit such large oscillations because a superposition of LCP and RCP 
light of substantially the same energy densities is needed to create them [33]. The LCP component of the 
refiected pulse is of short duration, on the order of the duration of the incident pulse. The pulse-bleeding 
effect only operates on the fraction of incident light in the Bragg regime that is RCP, and so the tail of 
the primary refiected pulse contains mostly RCP light but little to no LCP light when the circular Bragg 
phenomenon occurs. 

Corresponding plots for the nonlinear Cases Nl and N2 are shown in Figures [7| and |H1 respectively. 
While still small compared to the incident pulses', the peak values of U {z,t) in the heads of the primary 
refiected pulses increase from Figures IHl to |H1 for given carrier wavelength and polarization state. This is due 
to the increasing effects of nonlinearity in Case Nl as compared to Case L, and in Case N2 as compared 
to Case Nl. An understanding emerges from the frequency domain [10,15,33] as follows: for the type of 
nonlinearity we studied, a greater value of Ut creates a greater refractive-index mismatch between vacuum 
and the chiral STF because of the multiplicative character of the nonlinearity in (jSJ [33]. The refractive 
index of vacuum is unity, and we can see from Figure [21 that the real parts of n2 7^ 1 and 7^ 1 are 
decreasing functions of carrier wavelength at carrier wavelengths greater than 300 nm. When the chiral 
STF is nonlinear, the differences between the refractive indexes of the chiral STF and unity shall become 
even greater, and those differences will be greatest when the intensity of the beam is highest. Hence, the 
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reflection at the pulse head is greater than otherwise. 

The peak values of U (z, t) in the secondary pulses are not appreciably difl'erent from Case L to Case Nl 
to Case N2 (at least when normalized by the value of Ut, as was done for all the figures). Evidently, 
absorption by the chiral STF reduces the peak energy density of the primary refracted pulse in Case N2 
to a greater degree than in Case Nl and in Case Nl to a greater degree than in Case L [33]. In the 
frequency domain, absorption manifests itself through positive values of Im(n2) and Im(nrf), as seen in 
Figure IS By the time the primary refracted pulses reach zr, their energy densities drop such that the 
refractive-index mismatch between vacuum and chiral STF across z = is reduced between Cases L, 
Nl, and N2. Differences in energy density that remained were further reduced on the return trip of the 
secondary pulse from z = zr to z = z^. 

Figures lUl to ITTI contain plots of U {z,t) as functions of time recorded at z = 36 ^m; these are the 
transmitted pulses for Cases L, Nl, and N2, respectively. The peak values of U {z,t) increase, for given 
carrier polarization state, as Xcar increases, an effect that is attributable to the presence of the absorption 
resonances at ultraviolet wavelengths. The frequency domain behavior of the absorption resonances is 
shown in plots of Im (722) and Im (n^) in Figure [2 Other parameters being equal, more of the bandwidth 
of a pulse with lower Xcar overlaps the bandwidth of the absorption resonances of the chiral STF than the 
bandwidth of a pulse with higher Xcar, and hence more absorption occurs in the former case. For given Xcar 
and carrier polarization state, the peak values of U {z,t) — again normalized to Ut — drop from Figures 
to llll Notice that when the incident carrier plane wave is RCP, but not LCP, the transmitted pulses tend 
to have a lobed or humped structure, and that the number and sizes of the lobes increase from Case L to 
Case Nl to Case N2. Thus, the presence of the chosen nonlinearity distorts the shapes of the transmitted 
pulses and tends to make them highly asymmetric, but the effect is prominent only when the structural 
handedness of the film matches the sense of the circular polarization state of the incident carrier plane 
wave. 

3.3 Durations of transmitted pulses 

What effects does a change in carrier wavelength have on the durations of pulses transmitted through 
a chiral STF? Quantification of the effects of the carrier wavelength, carrier polarization state, and the 
nonlinearity of the chiral STF on Tp, t„, and Tc is provided through Figure I12( wherein plots of those 
quantities for thirteen values of Xcar ranging from 395 to 635 nm are shown. The incident carrier plane 
wave is again either LCP or RCP, and Cases L, Nl, and N2 are again considered. The overall trend in 
the plots is a decrease in pulse duration with increasing carrier wavelength. As the carrier wavelength 
decreases, the material exhibits greater dispersion due to proximity to the absorption resonances. The 
transmitted pulse is consequently of longer duration, because light of longer free-space wavelength has 
both greater phase and group velocities in the chiral STF than light of shorter Aq, as can be seen in 
Figure 01 This is because n2 (Aq) and rid i^o) drop with increasing Aq above the resonance wavelengths 
^^1,2,3; as can be seen in Figure |21 Therefore the chosen pulses, which do not possess any negative chirp in 
their carrier wavelengths, tend to spread out as they propagate through the chiral STF. 

Although that is the overall trend, there is an exception: there appears to be a local maximum for r„ 
when Xcar ~ 515 nm and the incident carrier plane wave is RCP. We believe that this effect is due to the 
circular Bragg phenomenon, because the local maximum appears in the vicinity of the center wavelength 
of the Bragg regime, and it does not appear when the incident carrier plane wave is LCP. Additionally, 
for the same Xcar and carrier polarization state, the duration of the pulses tends to increase from Case L 
to Case Nl to Case N2. 

3.4 Average transmission speed through chiral STF 

Plots of Cp, Cu, and Cc in FigureElshow that they all tend to increase with increasing carrier wavelength, but 
decrease slightly with increasing nonlinearity. These two trends are due to the presence of the absorption 
resonances in the ultraviolet regime, and the tendency of the nonlinearities to enhance those resonances. 
As is clear from Figure [2 the principal refractive indexes of the chiral STF are higher for wavelengths 
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closer to the resonances, and hence pulsed plane waves with lower Xcar ai'e slowed more by the film than 
pulsed plane waves with higher Xcar, as can be seen in Figure El 

The presence of nonlinearity increases the refractive indexes because of its multiplicative nature in ©, 
and hence pulsed plane waves with higher Ut are slowed more than pulsed plane waves with lower Ut when 
propagating through a chiral STF. Both Cp and Cc are heavily affected by pulse shaping by the chiral STF 
which is manifested as the alteration of the peaks of the transmitted pulses. For example, notice the jumps 
in these quantities that occur between Xcar = 455 nm and Xcar = 475 nm for Cases Nl and N2 when the 
incident carrier plane wave is RCP. The jumps occur due to the appearance of the lobed structure of the 
transmitted pulses mentioned in Section [3 .21 and, in the case of Cc, also due to the the short duration of 
the incident pulse as compared to the transmitted one. The calculations of Cp and Cc require finding the 
peak value of either U (zr, t) or of its correlation with the U (0, t); and so, if one lobe suddenly becomes 
larger than another, the peak of the transmitted pulse has to shift. This shift creates a jump in both Cp 
and Cc- Pulse shaping is to be expected to be particularly prominent when pulses of high peak powers pass 
through nonlinear chiral STFs, because the effects of nonlinearity are going to be very pronounced then. 

4 Concluding remarks 

As a result of our calculations, we offer the following conclusions that may be useful to future designers of 
optical pulse shapers fabricated with chiral STFs. 

• Increased nonlinearity reshapes pulses more — as is evident in the increased formation of a lobed 
structure — when the carrier polarization state of the incident pulse has the same sense as the struc- 
tural handedness of the chiral STF. Attenuation of the transmitted pulse increases when the carrier 
wavelength lies closer to the absorption resonances, which is also confirmed by frequency-domain re- 
sults [15,32]. 

• For specified carrier polarization state, Ut, and pnu the durations r^, r^, and Tc all tend to increase with 
decreasing Xcar — except for in the case of a linear chiral STF under conditions in which the circular 
Bragg phenomenon is exhibited. Therefore a chiral~STF~based device could change the durations of 
pulses, contingent on the polarization state of the incident carrier plane wave. 

• For specified Xcar and polarization state of the incident carrier plane wave, the transmitted pulse dura- 
tion tends to increase with increasing nonlinear effects. Therefore, any pulse-shaping device based on 
nonlinear chiral STFs can only be expected to work over a certain range of incident pulse intensities. 

• For a given carrier polarization state, Ut, and pni, the speed c„ tends to increase with increasing Xcar- 
Therefore chiral STFs can function as wavelength-selective delay elements for ultrashort optical pulses. 

• For specified Xcar and polarization state of the incident carrier plane wave, the speeds Cp, Cu, and Cc 
tend to decrease with increasing nonlinear effects. 

To conclude, let us note that our presented technique can accommodate other types of nonlinearities — 
specifically, nonlinearities that are not simply multiplicative but must be incorporated through multiple 
convolutions, which appear necessary for modelling the nonlinear response of dielectrics to attosecond 
pulses [20]. However, the computational capabilities then required will be significantly greater than the 
ones commonly available today. 
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Re(n2), Re(n^) 




Figure 2. Plots of Re (712) (top, dashed), Ke{n^) (top, solid), Im(n2) (bottom, dashed), and Im(n;j) (bottom, solid) as functions of 
free-space wavelength Aq (in nm). Note that these quantities decrease with increasing Aq for Aq > 300 nm. 





Figure 3. Plots of Vp2/cQ (top, dashed), v^j^/cq (top, solid), Vg^jcQ (bottom, dashed), and Vgii/cg (bottom, solid) as function of 
free-space wavelength Ao (in nm). Note that these quantities increase with increasing Aq. 
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Figure 4. Plots of U {z, t) (in arbitrary units) of the incident pulses as a function of z (in ^m) for \car = 415 nm (left), 515 nm 
(center), and 615 nm (right) recorded at time t = 63.4 fs. The incident carrier plane wave is either LCP (top) or RCP (bottom); note 
that the plots are identical regardless of carrier wavelength or polarization. The chiral STF, indicated by the shaded area, occupies the 
region 2 6 [30, 34] ^m. The values of U (2, t) have been multiplied by a factor of coToy/n/Ut; hence, the shapes of the incident pulses in 

these plots are independent of Ut also. 
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Figure 5. Same as Figure |4] except that the data are recorded at time t = 190 fs, and so reflected (left of film) and transmitted (right 
of film) pulses are captured. The chiral STF is structurally right— handed and linear (Case L). 
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Figure 6. Plots of U {z,t) (in arbitrary units) as a function of time t (in fs) for Xcar = 415 nm (loft), 515 nm (center), and 615 nm 
(right) recorded a,t z = 27.6 fim. The incident carrier plane wave is either LCP (top) or RCP (bottom). The primary and secondary 
reflected pulses are captured in the record. The chiral STF is structurally right-handed and linear (Case L). The values of U (z, t) are 

multiplied by a factor of coToy/n/Ut, where Ut = 1 x 10® J/m-^. 
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Figure 7. Same as Figure|S] except that the chiral STF is nonlinear (Case Nl). 
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Figure 8. Same as Figure[7| except that Ut = 2 X 10^ J/m^ (Case N2). 
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Figure 9. Plots of U {z,t) (in arbitrary units) as a function of time t (in fs) for Xcar = 415 nm (left), 515 nm (center), and 615 nm 
(right) recorded at 2 = 36 /im, and so the transmitted pulses are captured in the record. The incident carrier plane wave is either LCP 
(top) or RCP (bottom), and the chiral STF is structurally right-handed and linear (Case L). The values of U {z,t) are multiplied by a 

factor o{ CQToy/n/Ut, where Ut = 1 x 10® 
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Figure 10. Same as Figure|3 except that the chiral STF is nonlinear (Case Nl). 
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Figure 11. Same as FigureESl except that Ut = 2 X 10^ J/m^ (Case N2). 
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Figure 12. Durations Tp (A, in fs), (□, in fs), (x, in fs) as functions of carrier wavelength Xcar (in nm) for six different cases, 
evaluated at Zr = 36 ^tm over the interval [0, 190] fs. The incident carrier plane wave is either LCP (top) or RCP (bottom). In addition, 
the chiral STF is either linear (p„; 



b) while Ut = lx 10^ J/m^ (Case L, at left); nonlinear (p„i = 3 X 10-^4 m^/V^) while 



C/t = 1 X 10'' (Case Nl, in center); or nonUnear (p„; = 3 X lO""'** ra^ /Y^) while Ut = 2 X lO" J/m^ (Case N2, at right). The lines 

at the bottom of each plot indicate the value of Tp (dotted), Tu (solid), and Tc (dashed) for the incident pulses. 
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Figure 13. Normalized average speeds Cp/co (A), Cu/cq (□), and Cc/cq (x) as functions of carrier wavelength Xcar (in nm) for six 
different cases, evaluated at Zr = 36 ^m over the interval [0, 190] fs. The incident carrier plane wave is either LCP (top) or RCP 
(bottom). In addition, the chiral STF is either linear (p„; = 0) while Ut = 1 x 10^ J/m^ (Case L, at left); nonlinear 
(p„; = 3 X 10-2" m^/V^) while Ut = 1 x 10"^ J/m^ (Case Nl, in center); or nonlinear (p„; = 3 X lO'^" m^/V^) while 

Ut = 2x 10'^ J/m^ (Case N2, at right). 



